Electroweak Phase Transition in the MSSM: 4-Dimensional Lattice Simulations 
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Recent lattice results have shown that there is no Standard Model (SM) electroweak phase tran- 
sition (EWPT) for Higgs boson masses above ~ 72 GeV, which is below the present experimental 
limit. According to perturbation theory and 3-dimensional (3d) lattice simulations there could be an 
EWPT in the Minimal Supersymmetric Standard Model (MSSM) that is strong enough for baryoge- 
nesis up to m^RilOS GeV. In this letter we present the results of our large scale 4-dimensional (4d) 
lattice simulations for the MSSM EWPT. We carried out infinite volume and continuum limits and 
found a transition whose strength agrees well with perturbation theory, allowing MSSM electroweak 
baryogenesis at least up to mn = 103 ± 4 GeV. We determined the properties of the bubble wall 
that are important for a successful baryogenesis. 

PACS numbers: ll.lO.Wx, 11. 15. Ha, 12.60.Jv, 98.80. Cq 
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The visible Universe is made up of matter. This state- 
ment is mainly based on observations of the cosmic dif- 
fuse 7-ray background, which would be larger than the 
present limits if boundaries between "worlds" and "anti- 
worlds" existed 0|. The observed baryon asymmetry of 
the universe was eventually determined at the EWPT 
On the one hand this phase transition was the last 
instance when baryon asymmetry could have been gen- 
erated, around T w 100 - 200 GeV. On the other hand 
at these temperatures any B-l-L asymmetry could have 
been washed out. The possibility of baryogenesis at the 
EWPT is particularly attractive, since the underlying 
physics can be — and has already largely been — tested in 
collider experiments. 

The first detailed description of the EWPT in the SM 
was based on perturbative techniques , which resulted 
in 0(100%) corrections between different orders of the 
perturbative expansion for Higgs boson masses larger 
than about 60 GeV. The dimensionally reduced 3d ef- 
fective model (e.g. |^) was also studied perturbatively 
and gave similar conclusions. Large scale numerical sim- 
ulations both on 4d and 3d lattices were needed to ana- 
lyze the nature of the transition for realistic Higgs boson 
masses j5]j^. These results are in complete agreement, 
and predict |^,^ an end point for the first order EWPT 
at Higgs boson mass 72. Oil. 4 GeV ||^, above which only 
a rapid cross-over can be seen. The present experimen- 
tal lower limit of the SM Higgs boson mass is by several 
standard deviations larger than the end point value, thus 
any EWPT in the SM is excluded. This also means that 
the SM baryogenesis in the early Universe is ruled out. 

In order to explain the observed baryon asymmetry, 
extended versions of the SM are necessary. Clearly, 
the most attractive possibility is the MSSM. According 
to perturbative predictions the EWPT could be much 
stronger in the MSSM than in the SM [§, in particu- 
lar if the stop mass is smaller than the top mass pO[. 



At two-loop level stop-gluon graphs give a considerable 
strengthening of the EWPT (e.g. third and fourth paper 
of |]). A reduced 3d version of the MSSM has recently 
been studied on the lattice ^ (including SU(3) x SU(2) 
gauge fields, the right-handed stop and the "light" combi- 
nation of the Higgses). The results show that the EWPT 
can be strong enough, i.e. v/Tc>l, up to r7i/i«105 GeV 
and mj « 165 GeV (where rrih is the mass of the lightest 
neutral scalar and is that of the stop squark). The 
possibility of spontaneous CP violation for a successful 
baryogenesis is also addresed ||l^ . 

In this letter we study the EWPT in the MSSM on 
4d lattices and carry out infinite volume and contin- 
uum limit extrapolations. Except for the U(l) sector 
and scalars with small Yukawa couplings, the whole 
bosonic sector of the MSSM is kept: SU(3) and SU(2) 
gauge bosons, two Higgs doublets, left-handed and right- 
handed stops and sbottoms. As it has been done in the 
SM case fermions, owing to their heavy Matsubara 
modes, are included perturbatively in the final result. 
This work extends the 3d study pT| in two ways: 

a) We use 4d lattices instead of 3d. Note, that due to 
very soft modes — close to the end point in the SM — much 
more CPU time is needed in 4d than in 3d. However, this 
difficulty does not appear in the MSSM because the phase 
transition is strong and the dominant correlation lengths 
are not that large in units of T~^. Using unimproved 
lattice actions the leading corrections due to the finite 
lattice spacings are proportional to a in 3d and only to 
a? in 4d. For 0(a) improvement in the 3d case cf. [ p^ . 
In 4d simulations we also have direct control over zero 
temperature renormalization effects. 

b) We include both Higgs doublets, not only the light 
combination. According to standard baryogenesis sce- 
narios (see e.g. [Q) the generated baryon number is di- 
rectly proportional to the change of (3 through the bubble 
wall: A/3, (tan/3 = V2/V1, where vi^2 are the expectation 



1 



= 0.5 



0.4 



-1 — I — I — I — I — I — r- 



symmetric 



CB 



_L 



Higgs . 



J I I I I L 



I . < . I 



-1800 -1600 -1400 -1200 
m,= [GeV=] 

FIG. 1. The phase diagram of the bosonic theory obtained 
by lattice simulations. 

values of the two Higgses.) 

The continuum lagrangian of the above theory in stan- 
dard notation reads 



C — Cn + Ck + C\ 



Cy + Cw + 'Cs 



(1) 



The gauge part, Cg = 1/4 • F^I^Vt™)^"^ + 1/4 • FjfJ F'^^^^"' 
is the sum of weak and strong terms. The kinetic part is 
the sum of the covariant derivative terms of the two Higgs 
doublets {Hi,H2), the left-handed stop-sbottom doublet 
(Q), and the right-handed stop, sbottom singlets {U, D): 

{vj!^D*y{V'^''^f'D*). The potential term for the Higgs 
fields reads Cy = ml^[ai\Hi\^ + a2|i?2p - (^^1^2 + 

h.c.)] +gi/8 ■ + w - mi\w + mhIh,\^), 

for which two dimensionless mass terms are defined, 
Qfi = mf/m^2 s-iid ct2 = ^2/^12- One gets Csm — 
™qI*3P + "^cfl^l^ + "^dI^P foi' the squark mass part, 
and Cy = h^ti\QU\^ + |i?2p|C/|' + \Q^H2?) for the 
dominant Yukawa part. The quartic parts contain- 
ing the squark fields read Cw = ff^/8 ■ [2{(5}^ — 
\Q\^ + 4\hIq\^ + A\hIq\^ ~ 2|i7i|2|Q|2 _ 2|i72p|QP] and 
= 5sV8- [3{0r-|Q|4 + 2|C/|4 + 2|i?|4-6|Q[/|2 - 
6|QD|2 + 6|i7tD|2 + 2|Q|2|;7|2 +2\Q\^\D\'' - 2\U\''\D\''] , 
where {Q}^ = QiaQ^pQipQja- The scalar trihnear cou- 
plings have been omitted for simplicity. It is straightfor- 
ward to obtain the lattice action, for which we used the 
standard Wilson plaquette, hopping and site terms. 

The parameter space of the above Lagrangian is many- 
dimensional. We analyze the effect of the strong sector 
on the EWPT by using three specific sets of parame- 
ters. In one case the strong coupling has its physical 
value, whereas in the two other cases it is somewhat 
larger and smaller. The experimental values are taken 
for the weak and Yukawa couplings, and tan/3 = 6 is 



used. For the bare soft breaking masses our choice is 
iTiQ,D — 250 GeV, 771(7=0 GeV. Lattice renormalization 
effects on these masses will be discussed later. 

The simulation techniques are similar to those of the 
SU(2)-Higgs model Q (overrelaxation and heatbath al- 
gorithms are used for each scalar and gauge field); some 
new methods will be published elsewhere jl^] . The anal- 
ysis is based on finite temperature simulations (in which 
the temporal extension of the lattice Lt is much smaller 
than the spatial extensions Lx,y^z)i and zero tempera- 
ture ones (with Lt ~ L^^y^z)- For a given Lt, we fix 
all parameters of the Lagrangian except 0:2 • We tune 
Q!2 to the transition point, a2cj where we determine the 
jump of the Higgs field, the shape of the bubble wall, 
and the change of (3 through the phase boundary. Using 
a2c and the parameter set of the finite temperature case, 
we perform T — Q simulations and determine the masses 
(Higgses and W) and couplings (weak and strong) there. 
Extrapolations to the continuum limit and to infinite vol- 
umes are based on simulations at temporal extensions 
Lt — 2, 3, 4, 5 and at various lattice volumes for each Lt, 
respectively. Approaching the continuum limit, we move 
on an approximate line of constant physics (LCP), on 
which the renormalized quantities (masses and couplings) 
are almost constant, but the lattice spacing approaches 
zero. Our theory is bosonic, therefore the leading cor- 
rections due to finite lattice spacings are expected to be 
proportional to a? . This lattice spacing dependence is as- 
sumed for physical quantities in a — > extrapolations. 

We compare our simulation results with perturbation 
theory. We used one-loop perturbation theory without 
applying high temperature expansion (HTE). A specific 
feature was a careful treatment of finite renormalization 
effects, by taking into account all renormalization correc- 
tions and adjusting them to match the measured T = 
spectrum [|l^. We studied also the effect of the domi- 
nant r 7^ two- loop diagram ("setting-sun" stop-gluon 
graphs, cf. fifth ref. of [§), but only in the HTE. We ob- 
served less dramatic enhancement of the strength of the 
phase transition due to two- loop effects than in . Since 
the infrared behavior of the setting-sun graphs is not un- 
derstood, we use the one-loop technique with the T = 
scheme defined above. This type of one-loop perturba- 
tion theory is also applied to correct the measured data 
to some fixed LCP quantities, which are defined as the 
averages of results at different lattice spacings, (i.e. our 
reference point, for which the most important quantity 
is the lightest Higgs mass, ruh «45 GeV). 

Fig. |l]. shows the phase diagram in the m^-T plane. 
One identifies three phases. The phase on the left (large 
negative 777^^ and small stop mass) is the "color-breaking" 
(CB) phase. The phase in the upper right part is the 
"symmetric" phase, whereas the "Higgs" phase can be 
found in the lower right part. The line separating the 
symmetric and Higgs phases is obtained from Lt = 3 
simulations, whereas the lines between these phases and 



2 




the CB one are determined by keeping the lattice spac- 
ing fixed while increasing and decreasing the tempera- 
ture by changing Lt to 2 and 4, respectively. The shaded 
regions indicate the uncertainty in the critical tempera- 
tures. The phase transition to the CB phase is observed 
to be much stronger than that between the symmetric 
and Higgs phases. The qualitative features of this pic- 
ture are in complete agreement with perturbative and 
3d lattice results however, our choice of parame- 

ters does not correspond to a two-stage symmetric-Higgs 
phase transition. In this two-stage scenario there is a 
phase transition from the symmetric to the CB phase at 
some Ti and another phase transition occurs at T2 < Ti 
from the CB to the Higgs phase. It has been argued | |T^ ] 
that in the early universe no two-stage phase transition 
took place, therefore we do not study this possibiHty and 
the features of the CB phase any further. 

The bare squark mass parameters mQ,TO^,m|, re- 
ceive quadratic renormalization corrections. As it is well 
known, one-loop lattice perturbation theory is not suffi- 
cient to reliably determine these corrections, thus we use 
the following method. We first determine the position 
of the non-perturbative CB phase transitions in the bare 
quantities (e.g. the triple point or the T=0 transition for 
mfj in Fig. |l|). These quantities are compared with the 
prediction of the continuum perturbation theory, which 
gives the renormalized mass parameters on the lattice. 

Fig. H contains the continuum limit extrapolation for 
the normalized jump of the order parameter {v/Tc- up- 
per data) and the critical temperature (Tc/mw' lower 
data). The shaded regions are the perturbative predic- 
tions at our reference point (see above) in the contin- 
uum. Their widths reflect the uncertainty of our refer- 
ence point, which is dominated by the error of nih- Note 
that v/Tc is very sensitive to mh, which results in the 
large uncertainties. Results obtained on the lattice and 
in perturbation theory agree reasonably within the esti- 



mated uncertainties. (It might well be that the Lt—2 re- 
sults are not in the scaling region; leaving them out from 
the continuum extrapolation the agreement between the 
lattice and perturbative results is even better.) 

Based on this agreement we use one-loop perturbation 
theory without HTE to determine cosmologically allowed 
regions in the mj^ vs. nih plane of the full MSSM (in- 
cluding fermions, niA = 500 GeV), see Fig. |. The 
two lines for each m^^ (which intersect for lower val- 
ues of TTij^ ) correspond to upper bounds resulting from 
Vn/Tc — 1 (steeper curves, Bl) and the T=0 maximum 
MSSM Higgs mass (B2). w„ is the non-perturbative Higgs 
expectation value, assumed to be larger than the pertur- 
bative one by 14%, a correction factor obtained in the 
bosonic model (cf. Fig. ||). For large mg (e.g. 600 GeV, 
771(^=630 GeV) the region below Bl is below B2. De- 
creasing mg B2 decreases and Bl increases. At TOg=560 
GeV (mt^=590 GeV) Bl and B2 intersect at the CB 
value of mi^. Since B2 is almost constant this yields the 
overall maximum Higgs mass for a successful baryogene- 
sis. For even smaller mg (m^^) both Bl and B2 are rel- 
evant. Note that the maximum Higgs mass corresponds 
to a finite value of mg «560 GeV, yielding m/j ~ 103 ± 4 
GeV (including also the uncertainties from Fig. || and 
the difference between the one and two-loop maximum 
Higgs mass calculations |p7|). 

In order to produce the observed baryon asymmetry, 
a strong first order phase transition is not enough. Ac- 
cording to standard MSSM baryogenesis scenarios p^ ] 
the generated baryon asymmetry is directly proportional 
to the variation of (3 through the bubble wall separating 
the Higgs and symmetric phases. By using elongated lat- 
tice [2 ■ ■ 192), L=8,12,16 at the transition point we 
study the properties of the wall. In our simulation pro- 
cedure we restrict the length of one of the Higgs fields 
to a small interval between its values in the bulk phases. 
As a consequence, the system fluctuates around a con- 
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M28413/FKFP-0128/1997. The simulations were carried 
out on the 46G PC-farm at Eotvos University. 



FIG. 4. The profile of the bubble wall for both of the Higgs 
fields for the lattice 2 • L^i92. 



figuration with two bulk phases and two walls between 
them. In order to have the smallest possible free en- 
ergy, the wall is perpendicular to the long direction. We 
eliminate the effect of the remaining zero mode by shift- 
ing the wall of each configuration to some fixed position. 
Fig.^ gives the bubble wall profiles for both Higgs fields. 
The measured width of the waU is [A+B-\og{aLTc)]/Tc 
A=10.8±.l and B=2.1±.l. This behavior indicates that 
the bubble wall is rough and without a pinning force of 
finite size its width diverges very slowly (logarithmically) 
psf . For the same bosonic theory the perturbative ap- 
proach predicts (II. 2 ± 1.5) /Tc for the width. 

Transforming the data of Fig. || to |i?2 1^ as a function 
of l^ip, we obtain A/3 = 0.0061 ± 0.0003. The pertur- 
bative prediction at this point is 0.0046 ± O.OOIO. Thus 
perturbative studies such as [ p^ are confirmed by non- 
perturbative results. 

To summarize, we presented 4d lattice results on the 
EWPT in the MSSM. Our simulations were carried out 
in the bosonic sector of the MSSM. We found quite a 
good agreement between lattice results and our one-loop 
perturbative predictions. Using this agreement together 
with a careful analysis of its uncertainties, we determined 
the upper bound for the lightest Higgs mass for a success- 
ful baryogenesis in the full (bosonic-|-fermionic) MSSM, 
which turned out to be (103 ± 4 GeV) consistent with 
the 3d analysis (« 105 GeV). We analyzed the bubble 
wall profile separating the Higgs and symmetric phases. 
The width of the wall and the change in (3 is in fairly 
good agreement with perturbative predictions for typi- 
cal bubble sizes. Both the upper bound for nih and the 
smallness of A/3 indicate that experiments allow just a 
small window for MSSM baryogenesis. 

Details of the present analysis will be discussed in a 
forthcoming publication |T^ ]. 
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ian Science Foundation Grants OTKA-T22929-29803- 



[1] 
[2] 

[3] 

[4] 
[5] 

[6] 

[7] 

[8] 

[9] 



[10] 

[11] 
[12] 



[13] 
[14] 



A. G. Cohen, A. De Rujula, S. L. Glashow, Astrophys. 
J. 495 (1998) 539. 

V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, 
Phys. Lett. B155 (1985) 36. 

P. Arnold and O. Espinosa, Phys. Rev. D47, 3546 (1993), 
erratum ibid., D50, 6662 (1994); W. BuchmuUer et 
al., Ann. Phys. (NY) 234, 260 (1994); Z. Fodor and 
A. Hebecker, Nucl. Phys. B432, 127 (1994). 
K. Farakos et al., Nucl. Phys. B425, 67 (1994); 

A. Jakovac and A. Patkos, Nucl. Phys. B494, 54 (1997). 

B. Bunk et al., Nucl. Phys. B403, 453 (1993), Z. Fodor et 
al., Phys. Lett. B334, 405 (1994); Nucl. Phys. B439, 147 
(1995), F. Csikor et al., Nucl. Phys. B474, 421 (1996). 
K. Kajantie et al., Nucl. Phys. B407, 356 (1993); Nucl. 
Phys. B466, 189 (1996); O. Philipsen et al., Nucl. Phys. 
B469, 445 (1996). 

K. Kajantie et al., Phys. Rev. Lett. 77, 2887 (1996); 
F. Karsch et al., Nucl. Phys. B (Proc. Suppl.) 53, 623 
(1997); M. Giirtler et al., Phys. Rev. D56, 3888 (1997). 

F. Csikor et al., Phys. Rev. Lett. 82, 21 (1999) (for an- 
alytical end point results see second paper of H and 
W. Buchmiiller, O. Philipsen, Nucl. Phys. B443 47 
(1995)). 

G. F. Guidice, Phys. Rev. D45, 3177 (1992); J. R. Es- 
pinosa et al., Phys. Lett. B307, 106 (1993); A. Brignole 
et al., Phys. Lett. B324, 181 (1994); J. R. Espinosa, 
Nucl. Phys. B475, 273 (1996); B. de Carlos, J. R. Es- 
pinosa, Nucl. Phys. B503, 24 (1997); D. Bodeker et al., 
Nucl. Phys. B497, 387 (1997); J. M. Cline, G. D. Moore, 
Phys. Rev. Lett. 81, 3315 (1998) M. Losada, Nucl. Phys. 



B537, 3 (1999); |hep-ph/9905441 
M. Carena et al., Phys. Lett. B380, 81 (1996); Nucl. 
Phys. B524 3 (1998). 

M. Laine, K. Rummukainen, Phys. Rev. Lett. 80 5259 
(1998); Nucl. Phys. B535 423 (1998). 
K. Funakubo et al.. Prog. Theor. Phys. 99 1045 (1998); 
ibid., 102 389 (1999); M. Laine, K. Rumm ukainen, Nucl. 



Phys. B545 141 (1999); |hep-lat/9908045 . 
G. D. Moore, Nucl. Phys. B523, 569 (1998). 
P. Huet, A. E. Nelson, Phys. Rev. D53, 4578 (1996); 
Phys. Lett. B355, 229 (1995); M. Carena et al., Nucl. 
Phys. B503, 387 (1997); A. Riotto, Nucl. Phys. B518, 
339 (1998); Phys. Rev. D58 095009 (1998); N. Rius, 



V. Sanz, 


iep-ph/990746C; M. Brhlik, et al.. 


ph/9911243 





[15] 

[16] 
[17] 
[18] 
[19] 



F. Csikor, Z. Fodor, P. Hegediis, A. Jakovac, S. D. Katz, 
A. Piroth, in preparation. 

J. M. Cline et al., Phys. Rev. D60, 105035 (1999). 



M. Carena et al., hep-ph/0001002 



D. Jasnow, Rep. Prog. Phys. 47, 1059 (1984) 

J. M. Moreno et al., Nucl. Phys. B526, 489 (1998); 

P. John, Phys. Lett. B452, 221 (1999). 



4 



